Discrete and continuum spectra in the unified shell model approach 
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A new version of the nuclear shell model unifies the consideration of the discrete spectrum, 
where the results agree with the standard shell model, and continuum. The ingredients of the 
method are the non-Hermitian effective Hamiltonian, energy-dependent one-body and two-body 
decay amplitudes, and self-consistent treatment of thresholds. The results for helium and oxygen 
isotope chains well reproduce the data. 
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The standard nuclear shell model (SM) approach erects 
a wall between the description of intrinsic structure and 
reactions since it does not account for the continuum 
spectrum. This problem became acute due to experi- 
mental progress towards nuclei far from stability. The 
proximity of continuum in loosely bound nuclei influences 
all their properties which makes necessary a common de- 
scription of the discrete and continuum spectrum. The 
classical book formulated the unified approach to nu- 
clear reactions based on the SM. However, only recently 
realistic practical methods |M ^ 111 were suggested. 
Below we formulate a version of the shell model in con- 
tinuum (CSM) that naturally includes the conventional 
SM with discrete spectrum and gives in the same frame- 
work the description of reaction cross sections and decay 
channels with no restriction by single-nucleon decays. We 
show non-trivial applications to the chains of helium and 
oxygen isotopes. The formalism is based on the Fesh- 
bach method of projecting out the states |c; E) of 
"external" Q-spa.ce, corresponding to the decay channels 
c at continuum energy E, and constructing the effective 
Hamiltonian H that acts in the "internal" P-space of the 
many-body SM states 

Within the total space P + Q we solve the Schrodinger 
equation H\a) = E\a), where the full wave function \a) is 
a superposition of internal states |1) and external states 
\c;E). The elimination of external states leads to the 
closed equation for the internal part with the effective 
Hamiltonian 



n{E) = H + AiE)--WiE), 



(1) 



where the principal value term A is due to the off-shell 
processes of virtual excitation into channel space Q, 



(1|A(P)|2)^^P JdE 



, AliE'T A-^{E') 



E-E' 



(2) 



and the explicitly non-Hermitian term 



(1|W^(P)|2)=2^ AliEYA^^iE) (3) 

c (open) 



represents on-shell decays into the channels open at given 
energy. The amplitudes A'l{E) are the matrix elements 
(c; E\H\1) of the original Hamiltonian between Q and P 
spaces. The two new terms depend on running energy 
E so that we deal with a non-Hermitian and energy- 
dependent Hamiltonian. It is important that the observ- 
able quantities, such as the scattering matrix, can be 
written in terms of the same amplitudes and the Hamil- 
tonian 7i Q . The factorized form Q of follows from 
the unitarity of the S'-matrix. 

This formulation is exact. The complex poles of the an- 
alytical continuation of ^(i?), £a — Ea~{i/'2)Ta, are rcs- 
onances. Since the decay amplitudes Ai(E) must vanish 
below threshold energy their energy behavior com- 
plicates the situation making decay non-exponential. On 
the other hand, this ensures a continuous matching to the 
conventional SM: below thresholds the same bound states 
are obtained assuming that A{E) is included into renor- 
malization of the SM interaction. In this work we use 
m-scheme Slater determinants on bound single-particle 
orbitals in the mean field potential for P-space basis. 
The channel states are defined by their asymptotic quan- 
tum numbers. A one-particle channel of an A^-body sys- 
tem with quantum numbers j, energy of the particle e^, 
and the residual nucleus in the state \a;N ~ 1), so that 
E = Ea+€j, is labelled as |c) = 6](ej)|a; A^ - 1). The 
total spin-isospin quantum numbers are restored by the 
solution if the Hamiltonian respects this symmetry. 

We assume that one-body decays are determined by 
the single-particle part h of the full SM Hamiltonian, 

) = ^ a^, (e, ) (a ; - 1 1 1 1 ; , (4) 

where aj(e) describes the single-particle transition from 
the SM state v into continuum state j with energy e^, 
mediated by the mean-field interaction h. If the va- 
lence space is small, and each single-particle state is 
uniquely marked by spin-isospin quantum numbers, only 
one orbital v can couple to a given continuum channel, 
aj(e) — 5yja^{e); then v can be identified with a contin- 
uum index j. The non-Hermitian part of effective Hamil- 
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tonian becomes {E = £j + E^) 

{l\W{E)\2)=27r6,, ^ \aHe,r\{a; N - l\b,\l; N)\\ 

c (open) 

(5) 

Being in general a many-body operator, in some im- 
portant cases W is effectively reduced to a single-particle 
form. Thus, far from thresholds one can ignore the en- 
ergy dependence and use the closure to simplify the sum- 
mation in @, 

{l\W{E)\2)=27r6r2j2\^'\'\{l;N\b]b,\l;N)\\ (6) 
j 

As a result, W becomes a one-body operator that assigns 

I ■ I 2 

a width 7j = 27r a-' to each unstable single-particle 
state i — v and can be combined with the SM Hamil- 
tonian by introducing complex single-particle energies 
Ci/ = Cj/ — i^^/2. A similar picture emerges if the residual 
two-body interaction is weak, and the sum over daughter 
systems a in Eq. ((SJ is dominated by a single term. The 
single-particle interpretation of the continuum coupling 
W is generally valid when the removal of a particle docs 
not lead to a significant restructuring of the mean field. 
In the lowest order the phenomenon of decay does 
not change the structure of the internal wave function. 
The width of the many-body state a is then given by the 
expectation value Fq, = (a|H^|a); if is assumed to be a 
one-body operator, Fq — "fi/{ey){a\b\,hy\a). This defines a 
many-body decay width as a product of a single-particle 
width and the spectroscopic factor {a\b\,hy\(\) . 

The single-particle decay follows from the one-body 
scattering problem in an average potential V that en- 
ters the one-body Hamiltonian h; this determines the re- 
duced single-particle amplitudes Qi(e). Assuming spher- 
ical symmetry of V{r), the single-particle decay ampli- 
tude is found as a'(e) — Fi{r)V{r)ui{r) dr, where ui 
is the radial function, and (for a neutral particle) Fi (r) = 
{2fi/TTky/^krji{kr). For the near-threshold cases, the 
continuum admixtures are dominated by the long wave- 
length states; with the fc ^ behavior Fi{r) ~ (fcr)'+^, 
then a'(e) - e('+i)/2. 

The two-nucleon decay admixes the two-body contin- 
uum in the asymptotic form |c) — 6](e)6j,(e')|a; — 2). 
The channel state is characterized by energies of emitted 
particles, e and e'. With the one-body interaction used 
for the single-particle decay, the admixture from the two- 
particle continuum appears as a second order contribu- 
tion, 

^ . . 

that proceeds through intermediate states /? of a nucleus 
with A"— 1 particles. Eq. and corresponding contribu- 
tion to W, Eq. |j2Jl, allow for additional simplifications in 



the near-threshold region. The overall width behaves as 
7 ^ q'^+i+i ^ where q is the total available kinetic energy. 

In contrast to this sequential decay amplitude, a 
direct two-body transition requires the presence of a 
two-body interaction in the Hamiltonian. To describe 
this process, a pair amplitude is introduced, A1{E) ~ 
a^^\ei,e2) {a; N — 2\pl\1; N) where operator = {6^ (g) 
bi^'}L removes a pair, and only the quantum numbers L 
of the pair are conserved. In the long wavelength limit, 
the two-body decay amplitude scales with energy in ac- 
cordance with the phase space volume, identically to the 
sequential decay. The dominant contribution of orbital 
momentum L = results inW ^ q^; this "pairing" chan- 
nel is also favored by the short-range nature of residual 
forces. 

As a first application, we consider the chain of he- 
lium isotopes from ^He to ^'^He. The internal P-space 
contains two single-particle levels, P3/2 and pi/2- The 
interaction and single-particle energies are defined in 
d, 0- For the one-body channels, using the Woods- 
Saxon potential with the parameters adjusted for ^He, 
it was determined that, even for several MeV above 
threshold, the single-particle amplitudes can be described 
by the parameterization 73/2(e) = 0.608 e"^/^ MeV and 
7i/2(e) 0.3652 £3/2 MeV for the decay width from 
P3/2 and pi/2 states, respectively. The sequential two- 
body decay is computed using Eq. O with the near- 
threshold approximation for the energy dependence of 
single-particle amplitudes. The direct two-body decay is 
introduced only for the pair emission with total angu- 
lar momentum L = 0. It is assumed that all internal 
L = pairs couple to the continuum with the same am- 
plitude a(-^""^(ei, £2). The direct two-body amplitude is 
parameterized as a^-'""^^ (£1, £2) = (fi -I- £2)/3'y27r, where 
the numerical constant fixes the strength of the residual 
two-body interaction. This is the only parameter of the 
model, and it was adjusted to two-body decays of ^He. 
Figure n and Table HI show the results of the CSM calcu- 
lation in a good agreement with data. 

Several features are worth emphasizing, (i) By de- 
sign of the model with A{E) approximated by a con- 
stant and included in the adjusted SM Hamiltonian, en- 
ergies of bound states agree exactly with the results of 
the standard SM. The ground states of ''He, ^He and ^He 
are nucleon-stable in agreement with experiment, (ii) 
Energies of resonances deviate from the SM predictions 
— the continuum is restructuring internal states. For 
narrow separated resonances the effect is small. How- 
ever, in systems strongly coupled to continuum, internal 
phase transitions with formation of broad (super-radiant) 
and very narrow states can be found, see [lOj and refer- 
ences therein. A similar effect can be traced in Table |2 
so that decaying states are "pushed" further into con- 
tinuum, (iii) A detailed comparison reveals information 
about structure and dominant decay modes. The decay 
of the 2+ state in ^He is a sequential two-body process 
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FIG. 1: Results for He isotopes. For each isotope CSM re- 
sult is on the left, while experimental value is to the right; 
experimentally observed states are linked to theoretical coun- 
terparts. The decay width is shown on the right in the units 
of MeV. Energies (or centroids of decaying states) are not 
shown, but the energy scale is given on the vertical axis. 
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TABLE I; Comparison of the traditional SM and CSM with 
data for He isotopes. First two columns identify the mass 
number and spin of the state. Next five columns compare en- 
ergies: £(SM) traditional SM; E{a,) version of CSM with only 
one-body decays included; E{h) version of CSM with one- 
body decay and its second order contribution to the two-body 
process; S(CSM) full CSM including the direct two-body de- 
cay mode; i5(EX) experimental data, if available. Next five 
columns compare decay widths from CSM calculations with 
experimental values; the SM calculation gives only discrete 
energies. All numbers are given in MeV; energies are mea- 
sured from the ground state of '*He. 



leading to the ground state of '^He with no effects of di- 
rect L = 2 pair emission. This is not the case for the 0^ 
state, where one-body and sequential two-body decays 
reproduce only about half of the observed width. The 
rest comes from the direct two-body emission oi L ~ {) 
neutron pair to the ground state of ^He favored since 
the state is mainly a coherent excitation of a neutron 
pair to the pi/2 orbital. The largest deviation from ex- 
periment is seen for the 3/2~ state in ^He (although the 
data also have a large uncertainty). Possibly, the direct 



pair decay to 3/2~ state in ^He (as seen from Table, it is 
very important) is followed by a fast further breakdown 
to ^He. Sequential three-body decays were not consid- 
ered in the model. The admixture of the high- lying 3/2~ 
state with a large width is also possible, (iv) The results 
agree with calculations 0, by a different method but 
with a similar mean field and schematic residual force. 

As second application. Fig. |21 we perform a full cal- 
culation (all states and all interaction matrix elements 
in the sc?-SM included) for the chain of oxygen isotopes; 
the first results were reported earlier As for He iso- 
topes, self-consistency is established at all stages: decay 
energies for given parent-daughter pairs are consistent 
with the reaction processes; decay amplitudes for given 
states are in agreement with the reaction calculation and 
SM spectroscopic factors; the effect of decay on intrinsic 
states is accounted for in the diagonalization of the non- 
Hermitian energy-dependent Hamiltonian. The model is 
essentially parameter-free. The standard SM interaction 
(USD '11^ or its sHghtly modified version for heavier iso- 
topes HBUSD 12]) was used supplemented by a param- 
eterization of one-body continuum coupling found from 
scattering off the Woods-Saxon potential adjusted for 
^^O. Due to scarce experimental data on neutron decays 
in heavy isotopes, only sequential two-body processes are 
considered which requires no additional parameters. A 
comparison with data is shown in the inserts on Fig. El 
for the USD interaction. A number of features similar to 
those discussed for the helium isotopes can be noticed. 

We presented fully realistic CSM calculations and com- 
pared results with experiment. The main ingredients of 
the method are the SM with a good effective interac- 
tion in the discrete spectrum, the non-Hermitian effec- 
tive Hamiltonian for open channels, correct energy de- 
pendence of the decay amplitudes, and the self-consistent 
calculation of thresholds for a chain of isotopes. The high 
quality SM results for discrete states are reproduced by 
the CSM exactly. The comparison with experiment in- 
dicates similar high quality for resonances in the contin- 
uum. A model requires very few parameters, and most 
of them can be fixed via simplified calculations of scat- 
tering processes. The self-consistency between reactions 
and structure through energy dependence, thresholds, 
and parent-daughter relations is a part of the method. 
Consideration of continuum goes beyond single-particle 
processes. A detailed study of one- and two-body decays 
and comparison with data provide insight into structure 
of resonant states and interplay of intrinsic and contin- 
uum dynamics. 
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FIG. 2: CSM calculations for oxygen isotopes with the HBUSD interaction. States from yellow (long lifetime) to red (short 
lifetime) are resonance states. The insert on the upper right shows a more detailed picture for the lightest ^''O to ^^O isotopes. 
Decays from all states that are experimentally measured are shown with arrows. A full comparison between available data and 
the calculation is given for ^^O. Energies are expressed in units of keV. Comparison of widths with available data is given in 
the table in lower-left corner. For both inserts the interaction USD was used that works better in this mass region. 
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